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ABSTRACT 

The early-time optical spectrum of the Type Ic supernova (SN) 2002ap was 
characterized by unusually broad features, leading some authors to designate it a 
"hypernova." We present optical spectropolarimetry of this object 16 and 37 days 
after the estimated date of explosion. After correcting for interstellar polariza- 
tion, we find evidence for a high level of intrinsic continuum polarization at both 
epochs: p > 1.3% on day 16 and p > 1.0% on day 37. Prominent line polariza- 
tion is also seen, especially in the trough of the Ca II near-infrared triplet during 
the second epoch. When interpreted in terms of the oblate, electron-scattering 
model atmospheres of Hoflich (1991), our results imply an asymmetry of at least 
20% (day 16) and 17% (day 37). The data suggest a fair degree of axisymmetry, 
although the polarization angle of the dominant polarization axes are different 
by about 55° between the two epochs, implying a complex morphology for the 
thinning ejecta. In particular, there exists some spectropolarimetric evidence for 
a different distribution of Ca relative to iron-group elements. 

We also present flux spectra of SN 2002ap taken 131 and 140 days after the 
explosion. The spectra are characterized by a very weak continuum and broad 
emission lines, indicating that SN 2002ap has entered the nebular phase. The 
spectral features are substantially similar to those of "normal" SNe Ic in the 
nebular phase, and the emission lines are not significantly broader. However, 
some of the broad lines are very sharply peaked, and may possess a narrow 
component (probably unresolved by our spectra, FWHM < 400 km s _1 ) that 
is redshifted by ~ 580 km s~ x with respect to the systemic velocity of the host 
galaxy. 

Subject headings: polarization — supernovae: individual (SN 2002ap) — tech- 
niques: polarimetric 
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1. INTRODUCTION 

Although core-collapse supernovae (SNe) present a wide range of spectral and photo- 
metric properties, there is growing consensus that much of this variety is due to the state of 
the progenitor star's hydrogen and helium envelopes at the time of explosion. Those stars 
with massive, intact envelopes produce SNe of Type II-plateau (SNe II-P; II due to the 
presence of H in the spectrum, and "plateau" due to the characteristic shape of the light 
curve), those that have lost their entire hydrogen envelope (perhaps through stellar winds 
or mass transfer to a companion) result in SNe lb (Type I because of the absence of obvious 
hydrogen in their spectra), and those that have been stripped of both hydrogen and most 
(or all) of their helium produce SNe Ic; see Filippenko (1997) for a general review. 

Recently, a new subclass of objects has emerged whose members generically resemble 
SNe Ic (no hydrogen or obvious helium spectral features), but, unlike traditional SNe Ic, 
have spectra characterized by unusually broad features at early times, indicating velocities 
in excess of ~ 30,000 km s -1 . A few also possess inferred kinetic energies exceeding that 
of "normal" core-collapse SNe by more than a factor of 10 (see, e.g., Nomoto et al. 2001). 
These objects are colloquially referred to as "hypernovae," although not all of them are 
clearly more luminous or energetic than normal SNe Ic. There are currently 5 generally 
accepted members of this rare class: SN 1997dq, SN 1997ef, SN 1998bw, and, most recently, 
SN 2002ap and SN 2002bl (see, e.g., Iwamoto et al. 2000; Matheson et al. 2000b; Filippenko, 
Leonard, & Moran 2002; Mazzali et al. 2002; Kawabata et al. 2002). A related subclass 
of SNe exhibits many of the characteristics of these objects, but with hydrogen present in 
the spectra; the clearest examples are SN 1997cy and SN 1999E (see Filippenko 2001, and 
references therein), and they, too, are sometimes called hypernovae. 

Intense interest in hypernovae has been sparked not only by their peculiar spectral 
features, but also by the strong spatial and temporal association between the brightest and 
most energetic of these events, SN 1998bw, and the 7-ray burst (GRB) 980425 (e.g., Galama 
et al. 1998). This potential association has fueled the proposition that some (or, perhaps 
all) core-collapse SNe explode due to the action of a "bipolar" jet of material (Wheeler, 
Meier, & Wilson 2002; Khokhlov et al. 1999; MacFadyen & Woosley 1999), as opposed to 
the conventional neutrino-driven mechanism (Colgate & White 1966; Burrows et al. 2000, 
and references therein). Under this paradigm, a GRB is only produced by those few events 
in which the progenitor has lost all of its outer envelope material (i.e., it is a "bare core" 
collapsing), and is only observed if the jet is closely aligned with our line-of-sight (1-o-s). With 
such an explosion mechanism, one can well imagine that severe distortions from spherical 
symmetry may exist in the ejecta. 

Since all hot, young SN atmospheres are dominated by electron scattering (Wagoner 
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1981), which is an inherently polarizing process, spectropolarimetry provides a powerful 
probe of young SN morphology. The basic question is whether the SN is round: by symme- 
try, if the atmosphere is spherical the directional components of the electric vectors cancel to 
produce zero net polarization. In contrast, an asymmetric atmosphere will yield a nonzero 
polarization due to incomplete cancellation. From recent polarimetry of core-collapse events, 
a tentative trend has been identified: the degree of polarization (and, hence, asymmetry) in- 
creases with decreasing progenitor envelope mass (Wheeler 2000). Indeed, while the intrinsic 
polarizations of "normal" SNe II-P are generally below 1% (Leonard et al. 2001; Leonard & 
Filippenko 2001), a polarization of over 4% has been reported for an SN Ic, suggesting an 
axial asymmetry of more than 3:1 in this event (Wang et al. 2001). The number of events 
investigated in detail with spectropolarimetry, however, remains small. 

An additional technique that has been used to infer explosion asymmetry is the analysis 
of nebular line profiles. Recent simulations involving explosive nucleosynthesis in aspherical, 
jet-induced SN explosions predict that intermediate-mass and heavy elements such as iron 
are ejected (at high velocity) primarily along the poles whereas elements synthesized in 
the progenitor (e.g., He, C, Ca, O) are preferentially located at lower velocities near the 
equatorial plane in the expanding ejecta (Maeda et al. 2002; Khokhlov & Hoflich 2001; 
Hoflich et al. 2001). Maeda et al. (2002) model the effect that explosion asymmetry has on 
the spectra of hypernova explosions, and compare the results to the observed nebular line 
profiles of SN 1998bw. The main observable diagnostic for explosion asymmetry is found to 
be the ratio of the width of a probable [Fe II] blend near 5200 A to that of [O I] AA6300, 
6364 A. Maeda et al. (2002) find models of spherical explosions to be incapable of yielding 
velocity ratios greater than ~ 3:2, whereas aspherical models can generate ratios of more 
than 2:1. Although a specific degree of asphericity is not given, Maeda et al. (2002) find 
the nebular-phase line-width ratio of [Fe II] to [O I] for SN 1998bw to be inconsistent with 
a spherical explosion. 

In this paper, we report on spectropolarimetric observations obtained for the recent 
peculiar 1 Type Ic SN 2002ap in M74 (Kinugasa et al. 2002; Filippenko & Chornock 2002; 
Mazzali et al. 2002; Gal- Yam, Ofek, & Shemmer 2002) during the photospheric phase, and 
two nebular-phase optical flux spectra. We describe the observations in § 2, present our 
results and analysis in § 3, and discuss our conclusions in § 4. Additional details of the 
observations and analysis are given in Appendix A. As this paper was nearing completion, 
a submitted paper by Kawabata et al. (2002) describing similar spectropolarimetric data 



1 Bergcr, Kulkarni, & Chevalier (2002) recently used radio observations to conclude that SN 2002ap does 
not exhibit evidence for a large amount of relativistic ejecta, and hence is an "ordinary" SN Ic; however, we 
feel that it remains a somewhat peculiar SN Ic, given its unusually broad early-time spectral features. 
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appeared on astro-ph; we briefly compare our results with that work in § 4. 

2. Observations 

We obtained spectropolarimetry of SN 2002ap on 2002 February 14 and March 7 (16 
and 37 days, respectively, after the time of explosion derived by Mazzali et al. 2002), with 
the Low-Resolution Imaging Spectrometer (Oke et al. 1995) in polarimetry mode (LRISp; 
Cohen 1996 2 ) at the Cassegrain focus of the Keck-I 10-m telescope. After spending nearly 
three months in solar conjunction, SN 2002ap emerged from behind the Sun in June, and 
we obtained optical flux spectra of it on 2002 June 8 and June 17 (131 and 140 days, 
respectively, after the time of explosion) with the Kast double spectrograph (Miller & Stone 
1993) at the Cassegrain focus of the Shane 3-m telescope at Lick Observatory. On the 
second June epoch, spectropolarimetry of five distant Galactic stars near to the 1-o-s of 
SN 2002ap was also obtained in order to quantify the Galactic interstellar polarization (ISP) 
component contributing to the observed polarization of SN 2002ap. (Due to the very limited 
time of visibility, spectropolarimetry of SN 2002ap itself was not obtained on this date.) 
The slit position angle was set to within 3° of the parallactic angle (Filippenko 1982) for all 
observations in order to minimize differential light loss, and the total-flux spectra from all 
epochs were corrected for continuum atmospheric extinction and telluric absorption bands 
(e.g., Wade & Home 1988; Bessell 1999). A journal of observations is given in Table 1. 

The polarimetry data were reduced according to the methods outlined by Miller, Robin- 
son, & Goodrich (1988) and detailed by Leonard et al. (2001) and Leonard & Filippenko 
(2001). Additional details of the observations and analysis, including an investigation of 
the potential impact of second-order light and instrumental polarization on our results, are 
given in Appendix A. The polarization angle (P.A.) offset between the half-wave plate and 
the sky coordinate system was determined for all polarimetry epochs (at both Keck and 
Lick) by observing the polarized standard star BD +59°389 and setting its V-band polar- 
ization position angle (i.e., 9y, the debiased, flux- weighted average of the polarization angle 
over the wavelength range 5050-5950 A; see Leonard et al. 2001) equal to 98.09°, the value 
cataloged by Schmidt, Elston, & Lupie (1992). To check the precision of the P.A. offset at 
Keck, we also observed the polarized standard star HD 127769 on both nights, and found its 
P.A. to agree to within 2.4° of the values measured by Clemens & Tapia (1990) and Hiltner 
(1956). Since the values individually reported by the two catalogs disagree by 1.3° in P.A. 
and 0.14% in polarization level for this standard star, we do not ascribe much significance to 



instrument manual available at http://www2.keck.hawaii.edu:3636/. 
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the disagreement found between our measurements and the cataloged values. It is possible 
that this star exhibits polarization variability. The observed null standard star, HD 57702 
(Mathewson & Ford 1970), was found to be null to within 0.08% on both nights at Keck. 
(See Appendix A for details of the small amount of instrumental polarization that may have 
existed in one of our settings during the February 14 epoch.) At Lick, the P.A. offset was 
checked by observing HD 204827; its P.A., measured to be 58.1°, differs by only 0.6° from 
the value cataloged by Schmidt et al. (1992). The null standards observed at Lick, BD 
+32°3739 and HD 212311, were both found to be null to within 0.11%. 

3. Results and Analysis 

3.1. Observed Polarization 

The results for SN 2002ap are shown in Figures 1 and 2. Before estimating and then 
removing the ISP (§ 3.2), we make a few general remarks about the observed polarization. 
The nature of the polarization of SN 2002ap changed significantly between our two obser- 
vations. On day 16, the observed polarization rises from ~ 0.6% near 4000 A to > 1.5% at 
10,000 A, with the general rise punctuated by some features that may be associated with 
strong features seen in the total-flux spectrum. The P.A. also changes by about 15° across 
the spectrum. The notable feature in the polarization data from day 37 is an enormous po- 
larization change associated with the P-Cygni profile of the Ca II near-infrared (IR) triplet. 
The polarization changes between the two observational epochs, as well as the strong mod- 
ulations across spectral lines, indicate that some of the observed polarization is intrinsic to 
SN 2002ap, and is not interstellar in nature. In order to progress with further interpretation, 
however, we must attempt to estimate and remove the ISP from our data. 

3.2. Removal of ISP 

From reddening considerations we do not expect a large ISP, since any dust that polarizes 
the SN light should also redden it, and SN 2002ap is not highly reddened: E{B-V) ym = 0.07 
mag (Schlegel, Finkbeiner, & Davis 1998) and E(B— V)mia = 0.01 mag (Klose, Guen- 
ther, & Woitas 2002), for a total estimated reddening of E(B— V) f=a 0.08 mag. Since 
the empirical Galactic limit for the maximum observed polarization for a given reddening is 
ISP / E(B— V) < 9.0% mag -1 (Serkowski, Mathewson, & Ford 1975; see, however, Leonard 
et al. 2002a for a possible violation of this inequality for dust in an external galaxy), this 
implies an upper limit to the total (Galactic plus host galaxy) ISP of only 0.72%. 
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Table 1. Journal of Observations of SN 2002ap 









P.A. C 


Par. P.A. d 


Resolution e 


Range f 


Air 


Seeing 11 


Exposure 


Day" 


DT Date 


Telescope' 5 


(deg) 


(deg) 


(A) 


(A) 


MassS 


(arcsec) 


(s) 


16.3 


2002 Feb 14.2 


Kip 


76 


77 


12 


3400-10100 


1.5-1.6 


2.0 


2 x 450 


16.4 


2002 Feb 14.3 


Kip 


76 


75-76 


12 


3400-10100 


3.1-3.6 


4.5 


500 + 600 


16.4 


2002 Feb 14.3 


Kip 


76 


74 


12 


3930-8830 


1.9-2.6 


2.5 


4 x 600 


37.3 


2002 Mar 7.2 


Kip 


76 


74-75 


12 


3900-8800 


2.4-3.0 


1.5 


4 x 300 


37.4 


2002 Mar 7.3 


Kip 


76 


7.3-74 


12 


3900-8800 


3.3-4.1 


2.0 


200 + 200 + 240 + 300 


130.6 


2002 Jun 8.5 


L.3 


125 


125 


10 


3200-5500 


3.5 


7.0 


1200 


130.6 


2002 Jun 8.5 


L.3 


125 


125 


15 


5200-10400 


3.5 


7.0 


1200 


139.6 


2002 Jun 17.5 


L3p 


127 


124 


8 


4280-7040 


2.3 


4.0 


500 



Note. — The first two sets of Keck observations were performed in dual-beam mode with a D680 dichroic separating the blue and red 
parts of the spectrum; an OG570 order-blocking filter was also inserted on the red side. A 1'.'5 slit was used for all Keck observations, and 
a 3'.'0 slit was used for all Lick observations. Although the polarimeter was attached during the observation on June 17.5, only a total flux 
spectrum of SN 2002ap was obtained. 

a Day since estimated date of explosion, 2002 Jan 28.9 (HJD 2,452,303.4; Mazalli et al. 2002). 

b L3(p) = Lick 3-m/Kast Double Spectrograph (Miller &; Stone 1993; 'p' denotes polarimeter attached); Kip — Keck-I 10-m/Low 
Resolution Imaging Spectrometer (LRIS; Oke et al. 1995) with polarimeter. 

c Position angle of the spectrograph slit. 

d Parallactic angle range over the course of the exposures. 

e Approximate spectral resolution derived from night-sky lines. 

^Wavelength range of the calibrated flux spectrum. In some cases, the wavelength range for which the spectropolarimetry could be 
calibrated is somewhat less; see Appendix A for details. 
SAirmass range for each set of observations. 

h Average value of the full width at half maximum (FWHM) of the spatial profile for each set of observations, rounded to the nearest 075. 
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Fig. 1. — Polarization data for SN 2002ap obtained 2002 February 14, about 16 days after 
explosion (Mazzali et al. 2002). The NASA/IPAC Extragalactic Database (NED) recession 
velocity of 657 km s -1 for M74 has been removed in this and all figures, (a) Total flux, in 
units of 10 -15 ergs s _1 cm~ 2 A -1 , with prominent absorption features identified by Mazzali 
et al. (2002) indicated, (b) Observed degree of polarization, (c) Polarization angle in the 
plane of the sky. (d, e) The normalized q and u Stokes parameters, with the level of ISP 
estimated in § 3.2 indicated {dashed lines). The total-flux spectrum is shown at 2 A bin -1 , 
whereas the polarization data is binned at 20 A bin -1 to improve the signal-to- noise ratio. 
Unless otherwise indicated, the displayed polarization p is the "rotated Stokes parameter" 
(RSP; see Leonard et al. 2001) in all figures. 
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Fig. 2. — As in Figure 1, but for data obtained on 2002 March 7, about 37 days after 
explosion (Mazzali et al. 2002). Probable feature identifications are taken from Mazzali et 
al. (2002) and Patat et al. (2001). 
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Fig. 3. — Polarization data in the q-u plane for SN 2002ap on 2002 Feb. 14 {filled circles) 
and Mar. 7 {filled squares). Each point represents a bin 100 A wide, with the symbol size 
increasing from the blue end of the spectrum (A ~ 4000 A) to the red (10100 A for Feb. 
14, and 8800 A for Mar. 7). The best-fitting line to each epoch, determined by a uniform- 
weight least-squares fit, is drawn through the points; the point of intersection of the two lines 
represents a possible value of the ISP. When averaged with the results of another technique 
(see § 3.2), the adopted ISP ((frsp = — 0.17%,misp = —0.48%; open triangle) is derived. The 
general regions containing the points representing the polarization across the Ca II IR triplet 
are also indicated {absorption, emission). 
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One way to estimate the actual total ISP level is to rely on the data obtained for 
SN 2002ap itself, along with some theoretical assumptions. Figure 3 shows the observed 
polarization data in the q-u plane for both epochs. The majority of the data from each 
epoch appears to scatter about straight, but different, lines in the q-u plane. Such behavior 
is indicative of an axisymmetric asphericity in which the intrinsic P. A. is wavelength inde- 
pendent, a trait that has been identified in previous SN polarization studies (e.g., Howell et 
al. 2001; Leonard et al. 2001). If axisymmetry is suspected, then the ISP must lie along the 
best-fitting line to the points in the q-u plane; since the dominant axis for the two epochs 
prefer different P.A.s (~ 95° and 150° on days 16 and 37, respectively), the intersection of 
the two lines provides an estimate of the ISP. From Figure 3 we see that this occurs near 
giSP = -0.19%, wisp = -0.53%, or p lSP = 0.56%, ISP = 125°. Since ISP varies with 
wavelength according to a "Serkowski" curve (Whittet et al. 1992), it has a wavelength 
dependence; for our purposes, we take the derived ISP location to represent the value at 
Amax = 5500 A, the average wavelength of maximum ISP found from studies of Galactic 
stars. 

A second method to estimate the total ISP is to assume that some spectral region is 
intrinsically completely unpolarized, with all of the observed polarization therefore coming 
from ISP. There is both theoretical and empirical evidence that this may be the case in the 
Ca II IR emission profile on day 37. Since directional information is lost as photons are 
absorbed and reemitted in a line (Hoflich et al. 1996), broad, unblended emission features 
have traditionally been treated as intrinsically unpolarized in SN studies. In their spectropo- 
larimetric study of an early-time, emission-dominated spectrum of SN 1998S, Leonard et al. 
(2000) demonstrate that the basic paradigm of polarized continuum light being diluted by 
unpolarized emission lines is quite consistent with that data: the degree of depolarization 
across the broad-line features of SN 1998S changes in direct proportion to strength of the line 
at each wavelength across the profile. In the case of the Ca II profile of SN 2002ap, however, 
we see that this simple conceptual scheme fails, since the entire emission component, and 
even some of the absorption component, is polarized at roughly the same level (Fig. 2); such 
behavior is also seen in the Ha profiles of the Type II-P SN 1999em (Leonard et al. 2001). 
One possible interpretation of this is that not only are the emission-line photons unpolarized, 
but the underlying continuum photons are as well, having been absorbed and reemitted by 
the optically thick line. Such a situation is also described by Howell et al. (2001) for an 
SN la. If we accept this premise, and estimate the ISP from a 200 A wide bin centered on the 
emission peak at 8578 A, we derive gs578 = —0.12%, M8578 = —0.35%. Fitting a Serkowski 
law ISP curve (Whittet et al. 1992) to this value with an assumed wavelength of maximum 
polarization of A max = 5500 A, we derive <2te P = —0.15%, uj S p = —0.43% at A = 5500 A, 
or pisp = 0.46%, #isp = 125°. This result is quite similar to that derived earlier by the q-u 
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plane method, and is also consistent with the value Kawabata et al. (2002) found using a 
similar technique (pi SP = 0.64% ± 0.20%, ISP = 120° ± 10°). 

The consistency between our two independent total ISP estimates affords some confi- 
dence that we have determined an accurate value. Taking the average of the two results 
yields an estimate of the total ISP of gisp = —0.17%, wisp = —0.48%, or pi S p = 0.51% at 
#isp — 125°. Since the reddening from Milky Way (MW) dust is thought to dominate the 
total reddening of SN 2002ap, we would expect Galactic ISP to make up the bulk of the 
total ISP. To check this assertion, we obtained spectropolarimetry of five distant Galactic 
stars near to the 1-o-s of the object (the catalog of Heiles 2000 records no polarization mea- 
surements of any Galactic star within 1° of the 1-o-s to SN 2002ap, the nominal angular 
separation over which some consistency of ISP can be assumed) , and list the results in Table 
2. Although the values have some scatter, they are generally consistent with each other 
(they all fall in the same quadrant of the q-u plane) and confirm a non-negligible polariza- 
tion in this area of the sky from Galactic dust. Calculating the mean and la scatter of the 
individual values, we obtain qyrw = —0.11 ± 0.08%, mmw = —0.29 ± 0.14%, which yields 
Pmw — 0.32 ± 0.13% at #mw — 125° ± 12°. This is to be compared with the total ISP 
estimated by the earlier techniques of pisp = 0.51%, #isp = 125°. The general agreement of 
all of these values affirms our contention that Galactic ISP makes up the bulk of the total 
ISP, and also strengthens our belief that the total ISP estimated by the earlier techniques 
must be nearly correct. We therefore adopt pisp = 0.51% at 9i SP = 125° for A max = 5500 A, 
and proceed to remove this ISP from our data. 3 



3.3. Intrinsic Polarization 

The resulting ISP-corrected, intrinsic SN polarization and P.A. are displayed in Figure 4. 
The outstanding feature of the intrinsic polarization on day 16 is its strong increase with 
wavelength, from p « 0.2% near 4000 A to p « 1.3% at the cutoff near 10,000 A. Such a 
precipitous rise in polarization at optical and near-infrared wavelengths is completely unlike 
the Serkowski law dependence of ISP, and must be intrinsic to the SN. A similar redward 
rise in intrinsic polarization was seen in spectropolarimetry of the Type la SN 1999by, where 
it was attributed to the decreasing importance of line opacities, and the increased influence 
of continuum electron scattering at longer wavelengths (Howell et al. 2001). In particular, 



3 The shape of some of the Galactic stars' polarization spectra suggests that A max may be somewhat 
less than the nominal value of A max = 5500 A. Given the inherent difficulty of accurately measuring such 
low polarization values (see, e.g., Leonard et al. 2001), however, we choose to use the nominal value of 
Amax = 5500 A to correct our data. 



- 12 - 



Tabic 2. Probes of Galactic Dust Polarization Near the Line of Sight of SN 2002ap 







Distance 13 


Angular Separation 


IV 


Uy 


PV 


»v 


Star 


Spectral Type a 


(PC) 


(deg) 


(percent) 


(percent) 


(percent) 


(deg) 


BD +15°241 


F8 


> 170 


0.38 


-0.06 


-0.40 


0.40 


131 


BD +15°239 


F8 


> 160 


0.42 


-0.09 


-0.19 


0.21 


123 


BD +15°233 


FO 


> 210 


0.60 


-0.24 


-0.48 


0.54 


122 


BD +15°236 


F2 


> 200 


0.62 


-0.05 


-0.25 


0.25 


129 


HD 9510 


A2 


> 380 


0.74 


-0.12 


-0.14 


0.18 


115 



Note. — The typical statistical uncertainty in the measurement of qy and uy for these stars is < 0.05%; however, 
this does not incorporate the 0.1 — 0.2% systematic uncertainty that may exist in polarization measurements made with 
the Lick 3-m polarimeter; see Leonard et al. 2001. Qv,uy,py, and By represent the debiased, flux-weighted averages of 
the q Stokes parameter, u Stokes parameter, polarization, and polarization P.A., respectively, over the wavelength range 
5050 - 5950 A (see Leonard et al. 2001). 

a Taken from the Simbad Astronomical Database (http:simbad.u-strasbg.fr/Simbad). 

b Estimated from spectroscopic parallax, using the values from Binney & Merrifield (1998, p. 107) with the assumption 
that all stars are luminosity class V dwarfs. This assumption results in lower bounds on the inferred distances. In order to 
formally satisfy the criterion of Tran (1995) a MW star should be more than 150 pc above (or below) the Galactic plane 
to fully sample MW dust. The Galactic latitude of SN 2002ap (b ft! —46°) formally requires that stars be at least 209 pc 
away; two of our stars may be somewhat closer than this. 

c Between the probe star and SN 2002ap. 
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numerous lines, largely due to iron-group elements, dominate the spectrum at A < 6150 A 



and significantly depolarize the light, whereas the red end of the spectrum suffers less severe 
line-blending, allowing the intrinsically polarized continuum light to dominate. We suspect 
that a similar effect may be operating here, and that the high polarization (p m 1.3%) 
measured at the red end of the spectrum in fact represents a lower limit on the true level of 
continuum polarization undiluted by unpolarized line photons at this epoch. 

A detail of the region near the strong polarization increase associated with the absorption 
trough of the Ca II IR triplet on day 37 is shown in Figure 5. In the absorption trough there 
is an apparent change in observed polarization of Ap tot = 1.8%, where 



and Aq and Au are the individual changes in the normalized Stokes parameters compared 
with the nearby continuum value (note that Ap tot is an ISP-independent quantity; see 
Leonard & Filippenko 2001). It is possible to interpret this sharp polarization increase within 
the basic model first proposed by McCall (1984) and later extended by Jeffery (1991), that 
polarization peaks are naturally associated with absorption minima due to selective blocking 
of forward-scattered (and hence less polarized) light in P-Cygni troughs. With this simple 
model in mind, Leonard & Filippenko (2001) showed that a lower bound could be placed on 
the true intrinsic continuum polarization level of 



where Ap to t is defined by Equation 1, J cont equals the interpolated value of the continuum 
flux at the location of the line trough, and /trough is the total flux at the line's flux minimum. 
The main difficulty in the application of this to SN 2002ap is interpolating the continuum 
flux across the Ca II profile; our adopted estimate of the continuum level is indicated in 
Figure 5a. With Ap to t ~ 1-8% and J C ont /-f trough ~ 2.8, we find that p cont > 1.0% for this 
epoch. 

We therefore see evidence for significant intrinsic continuum polarization on both day 
16 and day 37 after explosion. When modeled in terms of the oblate, electron-scattering 
atmospheres of Hoflich (1991), the inferred continuum polarization demands an asphericity 
of at least 20% on day 16 and 17% on day 37; the actual value naturally depends on the 
(unknown) viewing orientation. It is also possible that ionization asymmetry or clumping in 
the ejecta may contribute to the polarization, which would reduce the inferred bulk asymme- 
try (see, e.g., Hoflich, Khokhlov, & Wang 2001). Detailed theoretical modeling of both the 



Ap tot = y/Aq 2 + Au 2 , 



(1) 



Pcont > 



Aptot 



(2) 
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Fig. 4. — Polarization level (top, thin lines) and P. A. (bottom) intrinsic to SN 2002ap after 
removal of an ISP characterized by pi^p = 0.51% at $isp = 125° for days 16 (left panels) 
and 37 (right panels) after explosion. For comparison of features, arbitrarily scaled total 
flux spectra (top panels, thick lines) are overplotted. Since the very low inferred polarization 
for much of the observed spectral range in the day 37 epoch renders the P.A. quite inde- 
terminate and an RSP difficult to define, we display the "optimal" polarization, defined as 
Ptrad - [(o- 2 (ptrad))/ptrad\, where p tra d = V V + « 2 (see discussion in Leonard et al. 2001, and 
references therein) for the day 37 data. 
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Fig. 5. — ISP-corrected spectropolarimetry of the region near the Ca II IR triplet for 
SN 2002ap on 2002 March 7, roughly 37 days after explosion. Error bars are la statis- 
tical for 20 A bin" 1 , (a) Total flux, normalized by a spline fit to the continuum, displayed 
at 2 A bin -1 for better resolution, (b) Intrinsic degree of polarization, (c, d) Normalized q 
and u Stokes parameters. The spurious flux change near 7600 A is due to the difficulty in 
accurately removing the A-band telluric feature from the high air-mass observations. 
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total flux and polarization data may set constraints on the viewing orientation (e.g., Hoflich, 
Wheeler, & Wang 1999), and also serve to check the specific assumptions that underlie our 
estimates of the intrinsic continuum polarization. 

3.4. Nebular Flux Spectrum 

Figure 6 shows the total flux spectrum of SN 2002ap obtained after the SN reappeared 
following solar conjunction; we have combined the observations taken 131 and 140 days after 
the explosion in order to improve the signal-to-noise ratio. The spectrum is characterized 
by strong emission profiles superimposed on a very weak continuum, implying that the SN 
has entered the nebular phase (Maeda et al. 2002; Sollerman et al. 2000). The greater 
strength of the lines relative to the continuum level in the spectrum of SN 2002ap compared 
with other SNe Ic at similar, although slightly earlier epochs (Fig. 6), is consistent with the 
rapid early spectral evolution reported by Mazzali et al. (2002), who found that SN 2002ap 
evolved at about twice the rate of SN 1997ef. 

Early-time spectra of SN 2002ap showed very broad spectral features, indicating high 
velocities in the outer ejecta (~ 30,000 km s -1 ; Mazzali et al. 2002). Given the extreme 
breadth of the early-time spectral features, one might expect the nebular-phase broad-line 
widths (taken as the FWHM of a Gaussian fit to the continuum-subtracted line profile; see 
Matheson et al. 2001) for SN 2002ap to be greater than those measured in other nebular 
SNe Ic. However, for [O I] AA6300, 6364 we measure 8600 km s -1 , consistent with the 
average value (7600 ± 2100 km s _1 ) of 9 SNe Ic reported by Matheson et al. (2001), and 
for [Ca II] AA7291, 7324, we measure 9800 km s _1 , also consistent with the Matheson et 
al. (2001) value (8700 ± 2700 km s _1 ). At this stage, there is evidently not significantly 
more line emission coming from high-velocity gas in this "hypernova" than is typical for 
"ordinary" SNe Ic. Similarly normal expansion velocities were also found at late times for 
the hypernova SN 1998bw (Patat et al. 2001). In addiion, the identification of spectral 
features extending out to nearly ~ 30, 000 km s" 1 in very early-time spectra of the Type 
II-P SN 1999gi (Leonard et al. 2002b) suggests that the presence of high-velocity material 
may be common among core-collapse SNe, but is only easily identified if very early-time 
spectra are obtained. 

Figure 7 shows the [O I] AA6300, 6364 and Mg I] A4571 line profiles for day 140 (shown 
instead of the day 131 spectrum due to its superior resolution). A striking feature of both 
profiles is the very sharp peaks to the flux distribution, which may indicate the presence of a 
narrow-line component, redshifted by about 580 km s _1 with respect to the systemic velocity 
of the host galaxy M74. Given the complex nature of the profiles, it is not certain whether 
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Fig. 6. — SN 2002ap (top) in the early nebular phase, compared with the peculiar Type Ic 
"hypernova" SN 1997ef (middle) and the "ordinary" Type Ic SN 19941 (bottom) at similar 
epochs with estimated day since explosion indicated. The spectra are scaled so that the 
strengths of the Ca II IR lines are approximately the same. Prominent emission lines are 
identified, and the zero- flux levels (dotted lines) are indicated. The [Fe II] blend near 5200 A 
may contain significant [O I] and [Mg I] emission at this early nebular epoch (see discussion 
by Maeda et al. 2002). Some of the narrow emission lines resulting from a superposed H II 
region have been clipped in the spectrum of SN 1997ef. 
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Fig. 7.— The [O I] AA6300, 6364 (thick line) and Mg I] A4571 (thin line) profiles of SN 2002ap, 
140 days after explosion. In addition to the broad lines, there are narrow emission spikes in 
both profiles, redshifted by ~ 580 km s^ 1 with respect to the systemic velocity of M74. 
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these narrow features are resolved in our spectrum, which has a resolution of Av 400 km 
s~ x at 6300 A. Patat et al. (2001) observed redshifted emission peaks in spectra of SN 1998bw 
at a similar phase. This emission could come from an excess of gas existing in lower-velocity 
regions close to the center of the SN nebula. It is also possible that it is coming from material 
previously lost by the progenitor star and now existing (in clumps?) in the circumstellar 
environment of the SN (e.g., SN 1999cq; Matheson et al. 2000a). Continued monitoring 
at sufficient spectral resolution is urged in order to trace the development of these narrow 
features. 

Given the spectropolarimetric evidence for asphericity during the photospheric phase 
(§ 3.3), it is worthwhile to investigate whether the nebular-phase line profiles show any 
evidence for asymmetry in the line-emitting regions at these later epochs. As discussed in 
§ 1, the main spectral diagnostic of asymmetry was found by Maeda et al. (2002) to be a 
large ratio (> 1.5) between the width of the [Fe II] blend near 5200 A and that observed in 
[O I] AA6300, 6364 (see Fig. 6 for the identification of these features). One difficulty with 
this diagnostic, however, is the complex nature of the [Fe II] blend's profile, which might 
include contributions from other lines (e.g., [O I] A5577; [Mg I] A5173) that can substantially 
widen the observed profile, especially during the early nebular phase. It is also difficult 
to estimate the "continuum" underlying this feature, since it evidently drops significantly 
across the profile. In their study, Maeda et al. (2002) analyzed the feature in a nebular 
spectrum of SN 1998bw taken 201 days after 5-band maximum, in which they felt confident 
that no extraneous lines contributed to the profile. Given the rapid spectral development 
of SN 2002ap, it is possible that this region in our spectrum, from ~ 130 days after B 
maximum, is similarly uncontaminated. Detailed spectral models are needed to confirm this 
assumption, but they are beyond the scope of this paper. If we assume that our feature is 
the same as that observed in SN 1998bw at the later epoch, then we derive a velocity ratio 
of ~ 2.8:1 for the features. This value is similar to that found for SN 1998bw, which would 
therefore imply an aspherical explosion for SN 2002ap in the same way that it was inferred 
for SN 1998bw. However, the relatively early phase of our spectrum, the difficulty inherent 
in measuring the width of the [Fe II] feature, and the potential for significant line blending 
make us view this preliminary result with some degree of skepticism. 



4. Discussion and Conclusions 

There is strong evidence from spectropolarimetry and tentative evidence from nebular- 
phase line profiles that SN 2002ap is not spherical. It is useful to compare the spectropo- 
larimetry of SN 2002ap with the small but growing number of core-collapse events for which 
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polarimetry measurements exist. The lower bounds on the continuum polarization level of 
SN 2002ap (1.3% and 1.0% on days 16 and 37, respectively) are larger than those typically 
derived for SNe II-P (Wheeler 2000), but not nearly as great as that inferred for the SN Ic 
1997X (Wang et al. 2001). 4 A more direct comparison to a similar event can be made with 
the spectropolarimetry presented by Patat et al. (2001) for the hypernova SN 1998bw. Patat 
et al. (2001) present spectropolarimetry from two epochs, —7 and +10 days after 5-band 
maximum (our observations of SN 2002ap occurred on days 7 and 28 relative to the date 
of B maximum derived by Gal- Yam et al. 2002) with coverage from 4000 A to 6900 A. 
Although the ISP is unknown for this event, 5 we note that the overall level of the observed 
polarization {p « 0.5 — 0.7%), and the marked increase in the polarization level from blue 
to red wavelengths (particularly in the earliest epoch), are quite similar to the behavior of 
SN 2002ap. 

In addition to the high inferred continuum polarization, another interesting finding 
for SN 2002ap is that while the data from both epochs suggest a substantial degree of 
axisymmetry (Fig. 3), the dominant axes have quite different polarization angles, changing 
from rsj 95° on day 16 to ~ 150° on day 37 (§ 3.2). This behavior is different from that 
observed in other core-collapse events, most notably those resulting from progenitors with 
massive envelopes intact at the time of explosion. For SN 1999em, a classic Type II-P 
event, Leonard et al. (2001) show that a common, wavelength-independent P.A. exists for 
all observational epochs, which includes days 12, 45, 54, 164, and 168 after explosion. This 
suggests a more complex morphology for the thinning ejecta of SN 2002ap. 

As mentioned in § 1, a potentially relevant result of recent simulations involving explosive 
nucleosynthesis in aspherical, jet-induced SN explosions is that intermediate-mass and heavy 
elements (e.g., Fe) as well as freshly synthesized radioactive material such as 56 Ni are ejected 
(at high velocity) primarily along the poles whereas elements synthesized in the progenitor 
(e.g., He, C, Ca, O) are preferentially located at lower velocities near the equatorial plane 
in the expanding ejecta (Maeda et al. 2002; Khokhlov k Hoflich 2001; Hoflich et al. 2001). 
Although the P.A. in the ISP-corrected data from day 37 (bottom, right panel of Fig. 4) 
is difficult to define in some regions due to the low polarization level, there is evidently a 



4 The presentation of data for SN 1997X by Wang et al. (2001) shows broadband averages of the po- 
larization data. It is therefore unknown what fraction of the claimed polarization (p > 4%) might be due 
to line-feature modulations similar to what is seen in the Ca II IR triplet of SN 2002ap during the second 
epoch, and not to intrinsically polarized continuum light. 

5 The Galactic ISP derived by adopting the measured linear polarization of a nearby star, HD 184100 
(Kay et al. 1998), is probably not appropriate since this star has a spectroscopic parallax of only ~ 14 pc, 
which is insufficient to fully sample the dust in the Galactic plane. 
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change near A 6000 A. Below this wavelength, 100°, a value similar to the one that 
characterizes the day 16 data. Above ~ 6000 A, 9 ~ 150° tends to dominate; this is especially 
easy to see (and believe, due to the high polarization) in the Ca II IR absorption trough 
(see also Fig. 3). Might this change in P. A. be due to the different spatial distributions of 
the iron-group elements, which dominate at blue wavelengths, from the elements synthesized 
by the progenitors, such as Ca and O, that are responsible for much of the opacity at red 
wavelengths (see Fig. 2a)? Along these same lines, an explanation for the nearly constant 
P. A. on day 16 at all wavelengths, and its similarity to the PA. of day 37 at blue wavelengths, 
could be due to the photosphere at early times being located in the highest- velocity (i.e., 
polar) material. By day 37, perhaps the SN has cooled sufficiently to expose the inner, 
slower-moving equatorial regions, and we are witnessing the competition between these two 
different distributions of material. 

To be sure, a small change in the ISP estimate, or even the wavelength of maximum 
polarization of the ISP, could alter the inferred P.A. of the continuum in the second epoch 
considerably. Further interpretation of the polarization change near 6000 A in the second 
epoch is probably not warranted. The ~ 55° P.A. difference between the the day 16 data 
and the region associated with the Ca II IR trough on day 37, however, remains regardless 
of the ISP choice, and may well indicate a true difference in the distribution of Ca relative 
to iron-group elements. 

Just prior to completion of this research, a paper describing similar data appeared on 
astro-ph by Kawabata et al. (2002). Although the observational epochs are slightly different 
(Kawabata et al. present data from 2002 February 9.2-11.3 and March 8.2-10.2), the general 
results are similar: an observed polarization level in the February epochs of p 0.6% that 
may rise somewhat toward the red (the increase is not as obvious as it is in our data since 
their data does not extend beyond 8300 A) and a dramatic polarization increase in the 
Ca II IR triplet's absorption trough followed by depolarization in the emission profile in the 
March data. One interesting speculation made by Kawabata et al. (2002) is that much 
of the polarized continuum at early times may result from scattering off of electrons in a 
jet, or bipolar jets, of material emitted from the SN during the explosion. Supporting this 
hypothesis is the general similarity demonstrated between the intrinsic polarized flux (i.e., 
^intrinsic x /) and the total flux spectrum redshifted by 0.23c, the presumed speed of the 
jet. We also find qualitative agreement between the intrinsic polarized flux and the total 
flux redshifted by z = 0.23c during our February epoch (Fig. 8), although we have not 
investigated how significant the correlation is. We do note, though, that the speculation by 
Kawabata et al. (2002) is consistent with the scenario proposed above that the polarization 
P.A. on day 16 is dominated by the distribution of the polar jet material. Countering the jet 
hypothesis, however, is the recent assertion by Berger et al. (2002) that the extreme faintness 
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Fig. 8. — The polarized flux of SN 2002ap on 2002 Feb. 14 (16 days after explosion; thin 
line) compared with the total flux spectrum from the same date redshifted by 0.23c (thick 
line), arbitrarily scaled and offset for comparison of features. 
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of SN 2002ap at radio wavelengths is inconsistent with the presence of a jet, regardless of the 
viewing geometry. In any event, additional spectropolarimetry of young SNe Ic is certainly 
warranted in order to further test the jet model. 

The total flux spectra obtained in the early-nebular phase confirm that SN 2002ap 
continues the rapid spectral evolution previously witnessed during the photospheric phase 
by Mazzali et al. (2002). It is characterized by a very weak continuum and broad emission 
lines, which demonstrate that it is already well along the transition to the fully nebular 
phase. The lines have widths similar to those observed in "normal" nebular SNe Ic. There 
are, however, unusual narrow lines superimposed on some of the broad-line profiles, including 
especially those of [O I] AA6300, 6364 and Mg I] A4571. It is possible that the ratio of the 
width of the [Fe II] blend near 5200 A compared with that of [O I] AA6300, 6364 is indicative 
of asymmetry in the line-emitting region, but additional, later-time spectra are needed to 
confirm this result. 
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Society's Small Research Grant Program. Our work was also funded by NASA grants GO- 
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by AURA, Inc., under NASA contract NAS 5-26555. Additional funding was provided to 
A.V.F. by NSF grant AST-9987438. We thank A. J. Barth and E. C. Moran for help with 
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suggestions that resulted in an improved manuscript. The assistance of the staffs at Lick 
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A. APPENDIX 

A.l. Details of the Observations 

The observational approach that we have taken at Keck in the past for our program to 
study SN geometry has been to obtain spectropolarimetry with a 600 groove mm -1 grating 
blazed at 5000 A, which yields a spectral range of about 4350-6850 A in the red arm of 
LRIS (Leonard et al. 2000; Leonard et al. 2001; Leonard & Filippenko 2001). The observing 
procedure of SN 2002ap on February 14 and March 7 differed in two ways from our traditional 
approach. First, in an effort to take advantage of the extended coverage possible with the 
recently commissioned blue arm of the LRIS double spectrograph (LRIS-B; an overview of 
the LRIS-B instrument is available at http:/ /alamoana.keck. hawaii.edu/inst/lris/lrisb. html), 
on February 14 we used the D680 dichroic filter to split the spectrum at 6800 A and send 
the blue and red light to the respective arms of LRIS (hereafter referred to as the "UV" 
and "IR" settings, respectively). To our knowledge, these are the first spectropolarimetric 
observations to be published using LRIS-B. Second, on both nights we observed SN 2002ap 
using the red side alone (i.e., no dichroic) with a 300 groove mm -1 grating blazed at 5000 A, 
which yielded coverage in the range 3930-8830 A (hereafter referred to as the "OPT" setting). 
The lack of an order-blocking filter in this setting makes second-order light contamination 
beyond ~ 6400 A a concern. In this Appendix, we discuss the performance of LRIS-B for 
spectropolarimetry and the impact of second-order light on the data acquired with the OPT 
setting. 

For the UV/IR setup, a grism with 300 groove mm" 1 blazed at 5000 A dispersed the 
light on the blue side, and a grating with 400 groove mm -1 blazed at 8500 A dispersed 
the light on the red side; an OG570 order-blocking filter was also inserted on the red side. 
In principle, such a setup yields coverage down to the atmospheric cutoff at ~3200 A. 
However, we found the calibration of the spectropolarimetry to be problematic at blue wave- 
lengths. At these wavelengths the poor efficiency of the blue CCD chip and grism (see 
http:/ /alamoana.keck. hawaii.edu/inst/lris/ccdblue_qe. html and 

http://alamoana.keck.hawaii.edu/inst/lris/grism_list.html), along with the inefficiency of the 
"HN-22 UV/optical Polaroid" polarizing filter, which is used to make the wavelength-dependent 
position angle correction curve (see Miller et al. 1988; Ogle et al. 1999), probably contribute 
to the difficulty we encountered. When we attempted to calibrate the polarization stan- 
dard BD +59°389 (Schmidt et al. 1992), strong, artificial PA. rotations were evident below 
~ 3800 A, indicating that the derived P. A. correction curve was not appropriate at these 
wavelengths. The data were also extremely noisy, especially for red objects. Fortunately, a 
new polarizer for the blue side that should allow calibration of polarization further into the 
blue has recently been installed and is currently being tested. In the end, we concluded that 
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reliable spectropolarimetry on the blue side was possible for SN 2002ap down to 3800 A. 
The total flux spectrum was calibratable down to 3400 A. 

An additional difficulty with the UV setting is the known problem of "ghost" images 
of second-order dispersed red-side light appearing on the blue side. Although the intensity 
of the "ghost" images is reported to be no more than 0.5% of the red-side intensity, such 
contamination, if variable between observations, could corrupt spectropolarimetry results 
in the UV setting. These uncertainties, as well as some previous problems that we have 
had calibrating spectropolarimetry data obtained with a dichroic (see, e.g., Leonard et al. 
2001), led us to approach the LRIS-B spectropolarimetry with some caution. The OPT 
setting fortunately provides a check on the UV data from this first epoch. In a similar way, 
the IR setting provides a check on possible second-order light contamination effects at red 
wavelengths in the OPT setting. 

On the red side, the "IR800 IR" polarizer, which is used to make the wavelength- 
dependent position angle correction curve, does not produce 100% polarized light at wave- 
lengths beyond ~ 8800 A. However, unlike the case on the blue side, consistent results (i.e., 
no unusual P.A. or polarization level changes) were obtained for all objects and standards out 
to 1 /an, indicating that the derived P.A. correction curve is accurate. This builds confidence 
in the calibration of the spectropolarimetry of SN 2002ap at red wavelengths as recorded by 
the IR setting. 

To check for instrumental polarization, the null standard HD 57702 (Mathewson & 
Ford 1970) was observed on February 14 with all three settings (Fig. 9). The polarization 
measured with both the OPT and IR settings are in good agreement and effectively null. 
Also impressive is the apparently good calibration of the IR polarization through the region 
affected by the dichroic cutoff. Similarly promising results were also obtained for the po- 
larization standard BD +59°389 (Schmidt et al. 1992) and other science objects observed 
on this night. Although the signal-to- noise ratio decreases sharply for both UV and IR 
near 6800 A, it may be possible to trust the polarization data in this region with no gap in 
coverage. 

The comparison between the UV and OPT settings for HD 57702, however, reveals that 
a small (< 0.1%) amount of instrumental polarization likely exists in the UV setting. This is 
also supported by similar discrepancies between the UV and OPT settings in other objects' 
data on this night. For HD 57702, Mathewson & Ford (1970) list p v = 0.04 ± 0.07%. From 
the OPT setting data we measure p v = 0.02 ± 0.02% (statistical error only), while the UV 
setting yields py = 0.10 ± 0.01%. Correcting for instrumental polarization at the 0.1% level 
is usually of dubious merit, since the intrinsic polarization of null standards is generally of 
this order. Nonetheless, since the OPT setting value is closer to the cataloged value, and we 
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Fig. 9. — Polarization of the null standard HD 57702 (Mathewson & Ford 1970) as recorded 
by the three settings described in the text, shown at 100 A bin -1 . Polarization degree was 
determined by the "optimal" polarization estimator (see caption to Fig. 4, and Leonard et 
al. 2001). 
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have never found significant instrumental polarization to be evident on the red side of LRIS, 
we attempted to correct for this small degree of instrumental polarization in the UV data. 
To do this, we first derived the offsets in the raw Stokes parameters between the UV and 
OPT measurements of HD 57702, and then used these offsets to correct the other objects' 
UV data (see Leonard et al. 2001). We found that this correction consistently improved the 
agreement between the two settings, although in some cases the agreement between the two 
settings was still not perfect. Since improvement was seen with the correction, though, we 
adopted it and applied it to the UV data obtained for SN 2002ap. The ultimate cause of 
this small amount of instrumental polarization in our UV setting is unknown. 

Although spectropolarimetry is routinely obtained at Keck using settings nearly iden- 
tical to our OPT setting (e.g., Brotherton et al. 2000; Ogle et al. 1999), this was the first 
time that we personally had obtained data significantly beyond 7000 A without the use of 
an order-blocking filter. This lead to some concern over the possibility of second-order light 
contaminating the data beyond ~ 6400 A. The degree to which second-order light affects 
red-side data has been investigated in detail by Putney & Cohen (1996) for the 300 groove 
mm -1 grating in the red camera, and they show that the percentage of blue light appearing 
as second-order red light in the range 6400-9000 A is at most 1.8% of the incident blue flux. 
Furthermore, from observations of polarized stars made with and without an order-blocking 
filter in place, Putney & Cohen (1996) conclude that second-order light effects have negligi- 
ble impact on polarimetry data, except perhaps for objects with highly wavelength- variable 
polarization (potentially a concern for SN 2002ap), and/or very blue flux. 

The redundancy of having the OPT data beyond 6400 A also covered by the IR setting, 
in which an order-blocking filter was used (in addition to the dichroic), allows us to directly 
test the conclusions of Putney & Cohen (1996) for SN 2002ap in the Feb. 14 data. Figure 10 
shows the polarimetry data for SN 2002ap obtained on this night in all three settings. The 
agreement at red wavelengths is indeed quite good, and suggests that second-order light 
effects are not a concern. At blue wavelengths, there is relatively good agreement between 
OPT and UV, although some discrepancies in the polarization level at the ~ 0.1% level do 
exist. 

Given the concerns with the UV data, the demonstrated lack of second-order light effects 
at red wavelengths in the OPT setting, and the very limited availability of SN 2002ap at 
the second epoch (it was approaching solar conjunction and visible for less than 1 hour 
immediately after sunset), we decided to observe SN 2002ap using only the OPT setting on 
March 7. Thus, direct comparison between these data and those for which an order-blocking 
filter was used is not possible for this epoch. However, we can estimate the fraction of the 
observed red flux contributed by second-order blue light through a procedure similar to that 
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Fig. 10. — Polarization of SN 2002ap obtained with the three settings described in the 
text, shown at 100 A bin^ 1 . UV setting data have been corrected for the small degree of 
instrumental polarization described in the text. 
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described by Putney & Cohen (1996), and compare this with the amount from the first 
epoch, for which second-order effects were shown to be negligible. 

Since the second-order light effects derived by Putney & Cohen (1996) were not derived 
using observations made with the polarimetry optics in the light-path, we determined our 
own values by observing the blue flux standard Feige 34 (Stone 1977; Massey & Gronwall 
1990) in the OPT setting, both with and without the GG495 order-blocking filter, as well 
as with the UV/IR setup. All observations were taken with the achromatic half-wave plate 
at a position angle of 0°, which produces two parallel, perpendicularly polarized beams on 
the CCD chip. The beams were individually flux calibrated and then combined to yield the 
total flux spectrum. Following the procedure of Putney & Cohen (1996), we calculate the 
fraction of blue light appearing as second-order contamination in the red end as 



h = (Al) 

^A/2 

where F' x is the flux-calibrated data measured from 6400-8828 A in the OPT setting without 
an order-blocking filter in place, F\ is the flux-calibrated data in the same region observed 
with an order-blocking filter, and F' x j 2 is the flux-calibrated data measured below 4414 A. 
We then smoothed f\ and interpolated over telluric features (e.g., the A-band) to produce 
the result shown in Figure 11, which is, in fact, quite similar to the result derived by Putney 
& Cohen (1996). 

To calculate the fraction of the observed red light contributed by second-order blue flux 
for SN 2002ap, then, we calculate 



TP F' x f 

-^secondorder A/2 J ^ {AO\ 

r observed r \ 

for both the February and March epochs, and show the results in Figure 12. (Flux-calibrated 
data were available for SN 2002ap down to 3400 A on February 14 and 3900 A on March 7; see 
Table 1.) It is clear that second-order light does not make an appreciably larger contribution 
to the March 7 data than it did on February 14. In both epoch's data, the contribution is less 
than 1.3% of the observed flux. This similarity, especially when coupled with the smooth, 
non- variable nature of the polarization in the range 3900-4400 A (Fig. 2), implies similarly 
negligible second-order light effects in the March data, at least at wavelengths beyond 7800 A, 
which includes the important Ca II IR region. 

Additional confidence that second-order light does not compromise the data from the 
March epoch is gained from the favorable comparison with the spectropolarimetry data that 



-30- 




Fig. 11. — Fraction of incident blue light appearing as second-order red light calculated using 
equation (Al) and data of Feige 34 taken with LRISp. 
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Kawabata et al. (2002) obtained at nearly the same epoch, presumably without any second- 
order light contamination. We therefore conclude that second-order light contamination has 
had a negligible impact on our data for both epochs, and concur with the conclusions of 
Putney & Cohen (1996) that second-order light is not a concern for the spectropolarimetry 
of objects that do not have either highly wavelength- variable polarization at blue wavelengths 
or extremely blue color. 
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Fig. 12. — Fraction of observed (total) flux contributed by second-order flux for SN 2002ap 
in the OPT setting observations on 2002 Feb. 14 and Mar. 7. 
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